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1.  Introduction 

The  problem  of  assigning  the  closed  loop  poles  of  a  linear  time 
invariant  multivariable  system  using  a  proper,  linear,  time  invariant, 
output  feedback  compensator  continues  to  be  of  great  interest.  Even 
though  several  outstanding  issues  remain,  good  progress  has  been  made 
as  evidenced  by  the  interesting  work  of  several  researchers.  Some 
recent  work  can  be  found  In'  the  references. 

In  particular  Kimura  1975,  Davison  and  Wang  1975  show  that  for  a 
controllable  observable  plant  (order  n,  m  outputs,  £  inputs)  it  is 
"almost  always"  possible  to  assign  min(n,m+£-l)  closed  loop  poles  arbi¬ 
trarily  close  to  a  given  set  of  real  and  complex  conjugate  values,  by 
using  constant  output  feedback.  The  issue  of  what  happens  to  the  remain¬ 
ing  unassigned  poles  is  not  addressed.  In  a  recent  paper  .Antsakl Is  and 
Wolovfch  1-977  present  a  different  way  of  assigning  m1n(n,£+m-l)  poles  and 
suggest  ways  of  dealing  with  the  remaining  unassigned  poles.  They  also 
extend  their  result  to  include  dynamic  output  feedback  and  show  that  if 
the  original  system  is  initially  augmented  by  q  Integrators  and  then 
constant  output  feedback  applied,  m+£+2q-l  poles  can  be  assigned. 

The  present  work  deals  with  the  question  of  how  many  poles  can  be 
assigned  when  the  order  of  the  compensator  is  fixed.  It  is  shown  that 
if  a  compensator  of  order  q  is  used  then  min(n+q,(q+l)m+q)  poles  can  be 
arbitrarily  assigned.  Throughout  the  paper  it  is  assumed  that  m>£.  This 


Is  not  a  restrictive  assumption  because  the  £2m  case  can  be  treated  In  a 
very  similar  way  and  "dual"  results  obtained,  (i.e.  min(n+q,(q+l)£+q  poles 

can  be  assigned).  The  method  of  attack  is  different  than  the  previous. .anas 
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result. 


It  is  further  suggested  how  the  locations  of  the  remaining 


unassigned  poles  can  be  controlled. 

The  work  proceeds  in  the  following  manner.  The  question  of  assign¬ 
ing  real  poles  is  addressed  initially.  The  Main  Lemma  takes  up  the  issue 
with  a  strictly  proper  plant  with  m  outputs  and  one  input.  This  result 
provides  insight  as  to  how  the  general  case  might  be  handled  and  is 
successfully  applied  to  the  mxt,  m£Jl  case  in  the  Theorem.  It  is  then 
shown  that  the  case  of  real  and  complex  conjugate  poles  can  be  treated 
in  the  same  manner. 
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2.  Formulation 

The  following  feedback  configuration  is  considered: 


u  I  _  v 


where  P  is  a  strictly  proper  mxa  input-output  transfer  function  which 
represents  the  plant  and  C  some  &)an  proper  dynamic  compensator.  Both 
P  and  C  have  elements  in  R(s)  the  field  of  rational  functions  in  s  over 
the  reals  R.  The  closed  loop  transfer  function  6  is  given  by 

6  -  P(I  +  CP)"1. 

If  the  following  notation  for  matrix  fraction  representations 
(Desoer  and  Vldyasagar  1975,  Kallath  198Q)  is  used: 

P  «  BRpARp  some  right  representation  of  P, 

*  ALPBLP  some  left  rePresentat‘ion  p» 

*  NRpDRp  some  right  coprime  representation  of  P, 

*  D[pNLP  some  left  coprime  representation  of  P, 

then  6  can  be  expressed  as: 

6  «  BRp(ALCARp  +  BLCBRp) 
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*  Wrp  +  nlcnrp^  \c  a  nrp  *  1dlc 


A#  A-1  A/ 

*  nrp  *  dlc 


(least  order). 


where  NRp,  4  are  right  coprime,  D^,  4  left  coprime. 

It  can  be  shown  (Chen  and  Hsu  1968)  that  if  $(s)  is  the  characteristic 
polynomial  of  the  closed  loop  system  then 


*(s)  =  a  det  4 

where  a  is  some  non-zero  constant. 


s1'1 
1  s 
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+  <t>„  we  know  that 
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r_1  _ii 

-  *52*  •  •  •  »*|(J 


where  Q  is  an  ixk  matrix. 

Several  definitions  of  genericlty  have  been  used.  Throughout  this 
paper  a  set  SCR*  will  be  called  generic  if  it  contains  a  non-empty 
Zariski  open  set  of  R*  (Zariski  and  Samuel  1958). 


-5- 


3.  Single-Input,  Multi pi e-Output  Case 

If  P  is  an  mxl  strictly  proper  transfer  function,  of  McMillan  degree  n, 
it  can  be  written  in  the  right  coprime  representation: 

|  P  *  N  d'1 

where 

K  d  =  s"  +  dn.j  s"'1  +  ...  +  dQ 

s"‘1  ♦  •••  +No  • 

Now  a  lxm  proper  compensator  of  order  q  can  be  written  in  the  left  represent¬ 
ation 

C  »  x_1  Y 

where 

|  *  -  »q  +  vi S<H  +  -  xo 

Y  «  Y  sq  +  Y,  ,  sq”1  +  ...  +  Y. 
q  q-1  0 

If  x,  Y  are  left  coprime  then  the  closed  loop  characteristic  poly¬ 
nomial  Is: 

4>(s)  *  xd  +  YN  ■  sn+q  +  ^n+q-1  s"^*1  +  ...  +  . 

This  relationship  can  be  expressed  In  the  following  way: 

Yq»  xq-i»  Yq-1»  •••»xo»  Yp^  ^q+l^d,N^  *  *  *  **o^ 


(3.1) 


where  Sq+1{d,N)  is  the  q+lth  order  Sylvester  Resultant  (Bttmead  et  al 
1978)  of  d»  N  Ca  (q-H ) (mxl )  x  (n+q+1)  matrix). 

From  the  above  one  can  clearly  see  that  the  coefficients  of  the  charac¬ 
teristic  polynomial  are  linear  functions  of  the  compensator  parameters. 

This  consideration  will  play  a  crucial  role  in  the  mxt  case  where  a  com¬ 
pensator  structure  will  be  employed,  that  satisfies  this  condition. 

Main  Lemma 

Let  P  *  Nd”1  be  an  mxl  strictly  proper  transfer  function,  with  d  of 
degree  n,  with  m|n  and  q>0  a  fixed  integer.  Let  k  *  (q+l)m  +  qs  n+q  *  i, 
and  4>(s)  be  the  closed  loop  characteristic  polynomial  and  define: 

W  »  {(N,d)  e  Rmn+n | d  -  sn  +  dn.1sn“1  +  ...  +  do,N  *  N^s""1  +...+  NQ} 


{(Sj,s2,  ...,  s^)  c  R  |s^  real} 

r/N  a  omn+n  v  For  which  there  exists  a  proper 
HN,a,sj  e  k  x  k  Qf  Qrder  such  that  s  s  • 

of  #(s)  12 


compensator 
sk  are  roots 


Then  Z  is  a  generic  subset  of  Rmn+n  x  Rk. 

Remark:  The  requirement  that  m|n  is  not  restrictive  and  it  is  introduced 
for  convenience.  The  general  case  can  be  treated  with  very  minor  alterations. 

Let  Sq+i»  be  the  submatrix  obtained  from  Sq+1(d,N)  by  removing  the  first 
row  and  first  column  and  a  the  first  row  of  Sq+1(d,N)  with  the  first  entry 
removed. 

Since  Sj,s2,.. . ,8^  will  be  the  roots  of  $(s)  iff 

Q  *  "  [Si1,s21,...,ski] 

S|,s2,..'.,S|(  will  be  the  roots  of  $(s)  if  we  can  find  a  £ 


such  that 


Z  *  Sq+l  *  Q  s  *  Cs^.Sg1 . ^3  -  a  •  Q  .  (3.2) 

The  compensator  which  will  accomplish  this  is 
C  «  x"1  Y  , 

x  *  sq  +  xq-1  sq_1  +  ...  +  xQ  , 

Y  =  Yq  sq  +  ...  +  Yq  . 

The  Main  Lemma  suggests  that  this  can  be  done  for  "almost  all"  N,  d  and 
s  *  (s,,  ....  sk). 

I 

Proof: 

It  is  required  to  show  that  Z  contains  a  non-empty  Zariski  open  set. 
The  matrix  Sq+^  Q  is  kxk.  Let  RrRmn+n  x  Rk  for  which  it  is  invertible. 

Vi 15  rank  k  on  a  Zariski  open  set,  which  is  non-empty  since  any  Nd"1 
with  equal  observability  indices  will  belong  to  this  set  CBitmead  1978). 

The  matrix  Q  is  also  rank  k  for  a  generic  subset  of  R  since  Q  is  rank  i: 
for  any  s,  for  which  the  s^  are  distinct.  The  product  will  be  invert¬ 
ible  for  a  generic  subset  of  Rmn+n  x  Rk.  Equation  (3.2)  therefore  does 
have  a  solution  for  a  generic  subset  of  Rmn+n  x  Rk. 

We  further  need  to  show  that  the  x(s),  Y(s)  so  constructed  are 
generically  left  coprime.  Clearly  x(s),  Y(s)  are  left  coprime  in  a  Zariski 
open'  set,  we  just  need  to  show  that  this  is  non-empty.  To  do  this  we  must 
suggest  an  (N,d,s)  which  can  be  thought  of  as  a  point  in  Rmn+n  x  Rk  space 


for  which  the  x(s),  Y(s)  are  left  coprime.  This  is  done  in  two  steps: 
a)  The  lxl  case. 

Let  S2.S2.-..»Sj(,S|(+1,...,Sj  be  real  and  distinct  with 

sl*s2*  •••  *si  >  0. 

Let  d(s)  *  sn. 

Let  ng,n2,...nn_2  be  defined  as  follows: 

(“1)  Og  *  ^2*^2*  ••• 

(-l)^”^  s  sum  of  all  possible  products  of  i-1  roots  at  a  time. 

1-2 

(-1)  ng  *  sun  of  all  possible  products  of  i-2  roots  at  a  time. 

(-l)*"n+*nn_1  =  sun  of  all  possible  products  of  i-(n-l)  roots  at  a  time. 
The  solution  of  (3.2)  then  becomes 

Z  ®  x^p  ...»  Xg,  Yg] 

where  Yq  «  Y  j  «  ...  *  Yj  *  0,  YQ  *  1 
and 

(-l)1_nXg  »  sum  of  all  possible  products  of  i- (n-l)-l  roots 
at  a  time 

» 

* 

2 

(-1)  xq_2  *  sum  possible  products  of  2  roots  at  a  time 

C-l)xq_i  ■  Sj  +  s2  +  s3  +  ...  +  s1  . 

An  s  Cs^ * Sg* - • -  »s1  >  0)  can  be  found  such  that  for  this  test 
point  Q  Is  full  rank  and  the  corresponding  x(s),  Y(s)  left 
coprime. 
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1 


b)  The  mxl  case. 

= 


.n-1 


Let  d(s)  =  s  +  s  +  ...  +  dQ 


t 


N(s) 


1 

1 

c 

1 — 

n-1 

Kn-2 

n-2  , 

1 - 

o 

1 _ 

_  nn-l_ 

s  + 

_  nr.-2__ 

s  +  . . .  + 

_no_ 

I 


1 


where  the  K-'s  are  m-lxl  vectors.  The  assignment  of 

J 

dp_jj...jdgj  n^_^,...,ng,  ^i’^2’****^i  is  as  in  part  a). 
The  Kj.'s  are'assigned  in  the  following  way: 

ol  q  n-q  n-q  n-1 

with  t  *  n-2q-l  let, 

Vi  “  *i 

Kq+2  =  K2 


Kq+t  *  Kt  ' 


Choose  Kj ,  <2 ,  ...  Kt  and  _s  i n  such  a  way  that 


Kt  Kt-1 


Kj  0 

k2  k2 


0  ...  Ka 


0 

0 


q+1  block  rows 


X  and  Sq+.|  Q  are  full  rank. 


The  solution  to  the  corresponding  equation  3.2  is: 
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Y' 


El 

i 

1 


S 


£  -  [yq>  xq.l . Xo’  U 

Yq  "  Yq-1  *  •••  ’  Y1  *  °-  Y0  *  CO,  0 . 1] 

The  xq_1>  ....  xQ  as  in  part  a). 

And  again  x(s),  Y( s)  are  left  coprime,  for  this  test  point. 

The  above  guarantee  that  x(s),  Y(s)  are  generically  left  coprime. 

This  completes  the  proof  of  the  Main  Lemma. 

The  Main  Lemma  suggests  that  for  "almost  all"  mxl  transfer  functions 
of  McMillan  degree  n  and  for  "almost  all"  s^  k  *  (q+l)m  +  q  there  exists 
a  proper  compensator  of  order  q  such  that  s^,s2»...,sk  are  k  roots  of  the 
closed  loop  characteristic  polynomial  <fr(s).  Since  <|>(s)  is  of  degree  n+q 
this  means  that  in  general  there  are  n-(q+l)m  unassigned  poles.  This  means 
that  some  of  these  could  even  be  unstable.  As  this  is  a  matter  of  concern 
it  is  presently  under  investigation.  It  is  possible  that  by  restricting 
the  assignable  roots  to  lie  in  a  certain  region  to  assure  that  all  of  them 
are  stable.  In  the  general  mxz  case  it  is  shown  how  additional  compensator 
parameters  can  be  introduced  to  help  control  the  remaining  unassigned  poles. 

In  the  last  section  it  will  be  shown  that  real  and  complex  conjugate 
values  can  also  be  considered  and  that  the  Main  Lenina  continues  to  hold. 
Remark:  It  is  interesting  to  note  how  the  number  of  assignable  poles 
increases  as  a  function  of  q,  the  order  of  the  compensator.  Assuming  that 
m|n  and  since  "generically"  the  observability  indices  of  the  plant  are  all 
equal  to  v  =  p  we  see  the  following: 

If  q*0,  m  poles  are  arbitrarily  assigned. 
q*l,  2m+l  poles  are  arbitrarily  assigned 

q*y-l,  iim  +  y-1  «  n+y-1  poles  are  arbitrarily  assigned, 


which  means  that  all  closed  loop  poles  can  be  arbitrarily  assigned.  A 
similar  result  when  q«y-l  has  been  obtained  earlier  by  Brasch  and 
Pearson  T970using  a  different  approach. 

An  algorithm  for  constructing  the  solution  involves  expressing  P  in  a 
right  coprime  representation  and  solving  a  linear  equation  over  the  reals 
(or  complexes)  since  £  from  (3.2) 

X  -  -  [sj,  . sJ](Sq+1  •  Q)"1  -  aQ(Sq+1Q)_1  • 

Comparing  this  procedure  with  the  one  suggested  by  Antsaklis  and 
Wolovich  1970  for  the  mxl  case  one  can  see  that  for  a  compensator  of  order 
q,  using  the  present  method  (m+l)q  +  m  poles  are  assigned,  whereas  with  the 
earlier  one  2q  +  m  are  assigned.  It  should  also  be  pointed  out  that  the 
method  used  there  is  different  than  the  present  one  in  this  respect  as  well 
in  that  initially  P  is  augmented  by  q  integrators  and  then  constant  output 
feedback  is  used  to  close  the  loop. 

The  example  below  helps  to  illustrate  and  clarify  the  procedure. 


a)  Let  q*0,  i.e.  a  constant  compensator.  The  result  suggests  that  2  poles 
can  be  arbitrarily  assigned. 

Let  Sj  *  -1,  s2  *  -2. 

The  compensator  used  is  given  by 


x  *  1  Y  »  [yx,  y2] 
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and  y^  y2  are  obtained  as  the  solution  of 


fvyzl  •  p  0  0  0  1  l 

011001 


r  5  Sr 

S1  s2 

4  4 

S1  s2 

3  3 

S1  s2 

.2.2 


-  -  Cs16,s263 


-  [0,  0,  0,  1.  0,  2]  .  Q 


Computing  the  solution  yields 

yl  =  16  ’  y2  =  '4 
and  the  compensator 

c-tU  .  -«]. 


b)  Let  q=l.  Then  5  poles  can  arbitrarily  be  assigned. 
Let  Sj  *  -1,  s2  *  -1.5,  s3  =  -2,  s4  *  -2.5,  s5  =  -3. 
The  compensator  used  is  given  by 

x  *  s  +  xx  ,  Y  *  [y1>  y2]s  +  [y3,  y43 


It  is  obtained  as  the  solution  of: 


» 


\ 

s 

1 


1 
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Cy i »y2»xi»y3»y4^ 


1  0  0  0  1  2  0 
0  110  0  10 
1  0  0  0  1  0  2 
0  1  0  0  0  1  2 
0  0  1  1  0  0  1 


1111 


Q 


*  -  [s17,s27,s37  s47  s57]  -  [0,  0,  0,  1,  0,  2,  0]Q  . 


Computing  the  solution  yields 
Xj  *  5.242544771594, 

y1  =  3.436523038651  y2  *  24.13305548336 
y3  =  2.405828042428  y3  *  7.162876396989 


and  the  compensator 

c »  pis  *  y3  • 

s  +  Xj  *  S  +  Xj 

c)  Let  q=2.  Then  all  8  poles  of  the  closed  loop  system  can  be  assigned. 
Let  Sj  *  -1,  s2  =  -1.1,  s3  =  -1.2,  s4  =  -1.3,  sg  *  -1.4,  sg  *  -1.5, 

Sy  *  -1.6,  Sg  =  -1.7. 

The  compensator  which  accomplishes  this  is  given  by  C  a  x-1Y 


x  «  sz  +  XjS  +  x2  Y  «  [ylt  y2]sz  +  [y3,  y4]s  +  [yg,  yg] 


where 


X1 

*1 

*3 


275.0363186229 

-335.7496784115 

-91,31887959451 


x2  «  63.27143983832 
y2  -  68.68243973376 
y4  -  -41.2992798809 
v.  «  17.61171990276 


v-  «  -6  766760056678 


4.  Multiple-Input,  Multiple-Output  Case 


If  P  is  an  mxt  (m>t)  strictly  proper  transfer  function  of  McMillan 
degree  n  and  equal  controllability  indices  X,  (i.e.  n*u)  it  can  be 
expressed  in  the  right  coprime  representation: 

p  *  Vrp 

where 

drp-  IsX  +  D*-i$k'1  *  tD0  <«■» 


— 1 
z 
>* 

1 

_ 1 

x-1 

"V 

3-i_ 

s  +  ...  + 

-K0. 

D.. ,  fL  are  fcxa  matrices,  K..  are  (m-t)xt  matrices,  flow  if  an  £xm 
proper  compensator  C  =  X”1  Y  is  used  of  the  form 

x  ■  V  *  v/'1  +  •••  +  xo 

Y  ■  V" +  Vi5"'1  +  -  +  Yo 

(X^  are  m,  Y^  are  txm),  and  X,  Y  are  left  coprime  then  the  closed  loop 
characteristic  polynomial  is  given  by 

♦(s)  *  det(X0Rp  +  YNRp)  . 

■  $ 

If  one  uses  Sylvester  Resultants  to  express  the  relationships 
XDRp  +  YNRp  one  has  the  following: 
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i 


*0>  U 


!*  B\-V  °X-2* 


0  ^x-1*  ^x-2* 


0  KX-1’  KX-2’ 


I 

0 

0 


•  •  •  »  Dg,  0  » 

"o’  °* 

•  •  • 

"o’  °’ 

DX-l’DX-2* 

^X-l*^X-2* 

KX-l*KX-2* 


,  0“ 
,  0 
,  0 

»  0 
,  0 
,  0 


q+1  block 
rows 

(4.2) 


*  £*X+q*  ^X+q-l . * 

where  4>  *  *x+(^x+<*  +  ...  +  *0>  is  fcxi. 

Now  since  <f>(s)  *  det  *,  it  is  evident  that  in  general  the  coefficients  of 
*(s)  are  non-linear  expressions  in  the  compensator  parameters.  This  is 
precisely  where  many  difficulties  concerning  the  pole  assignment  problem 
lie.  One  way  of  proceeding  is  to  find  ways  of  exploiting  this  non-linear 
structure.  In  this  paper  a  successful  approach  is  presneted  which  ’’avoids” 
the  non-linear  analysis.  The  problem  is  formulated  in  such  a  way  that  the 
nonlinear  structure  is  forced  to  become  linear. 

Suppose  that  the  compensators  under  consideration  are  restricted  to 


have  the  following  structure. 


Therefore 


x(s),  0,  ...  0 
0 

I.  , 


yu(s) »yi2(s)»  •••  yu(s),  ... 
0 


By  construction  C  *  X_1Y  Is  proper. 

Now  If  x(s)  is  the  first  row  of  Y,  (dn(s)  ,di2(s) . d^fs))  the  f1rst 

row  of  DRp,  £j(s)  the  column  of  NRp  and  (♦u(s),^12(s)>...*u(s)) 

the  first  row  of  ♦  one  has  the  following: 

x(s)du(s)  +^(s)ni(s)  ■  <^i  1  ( s ) 

x(s)di2(s)  +  ^.(s)ji2<s)  -  4> j 2 ( s ) 
x(s)du(s)  +J<(s)nt(s)  •  *u(s)  . 


Let  ?  be  the  (£-l)x£  matrix  which  contains  the  remaining  £-1  rows  of  #. 

Then 

•  =  [0,  I£.1]sX  +  (D'x_1+Nx_1)sX‘1  +  ...  +  (Dq-hTq) 

where  D\,  F.  are  obtained  from  »  N..  by  removing  the  first  row. 

Now  since  the  closed  loop  characteristic  polynomial  is  <|>(s)  *  dett  one 
can  easily  see  by  expanding  the  first  row  that: 

♦(s)  *  +  *12^^12  +  ***  +  *l£^Al£  (4*4) 

(A^  is  the  appropriate  (£-l)x(£-l)  minor  of  4>).  The  coefficients  of  <t>(s) 
are  linear  expressions  in  the  compensator  parameters.  This  relationship 
is  made  more  explicit  in  the  following: 

♦(s)  -  x(s)(d11(s)A11  +  ...  +  du(s)Au)  +  ^(s) (n.1(s)A11  +.. .+n£(s)Al£) 

*  x(s)d  +  ^(s)N.  (4.5) 

Since  d  is  a  polynomial  of  degree  X£*n,  N  an  mxl  vector  of  degree  A£-l*n-l, 
x(s)  a  polynomial  of  degree  q  and  y,(s)  a  lxm  vector  of  degree  q  this  fits 
precisely  the  formulation  used  in  the  analysis  of  mxl  systems.  Namely 
(4.5)  can  be  expressed  as: 

£l»Yq»xq-l*  •••»  xo,Y0^q+l^,N^  *  ^’*n+q*  ***  *0^  (4*6) 

This  would  Indicate  that  the  Main  Lemma  could  somehow  be  used  in  obtain¬ 
ing  a  result  for  the  mx£  case.  In  a  sense  the  original  mx£  system  has  been 
"transformed"  (or  reduced)  to  the  mxl  system  corresponding  to  (4.5).  This 
"transformation"  will  play  a  crucial  role  in  the  proof  of  the  upcoming 
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Theorem,  and  is  therefore  now  made  more  precise. 

Let  the  given  transfer  function  (4.1)  be  parameterized  as  follows: 


°RP  * 

- 

0 

• 

sX+ 

ai*  ***a£ 

aU+l’-**a(X+l)A 

• 

sx_1+...+ 

a(x-i)£+i**-ax£ 

a2X£ 

• 

0 

■ 

• 

1 

a(£-l)X£+l'‘* 

®£*X*£ 

'x-1 


NRP  = 


m-M 


d  A  •  •  •  3  A 

X*  +1  X£  +£ 

• 

dp  •••dp 

X£  +(X-1)£+1  X£  tX£ 

• 

• 

X£2+(£-l)X£+l 

sx‘*+...+ 

• 

• 

•  •  •  dp  p 

•  •  •  dp 

2X£  +X£ 

• 

• 

a  2 

2X£  +1 

• 

• 

• 

dp  •  •  • 

• 

•  •  •  dp 

_2X£  +(m-£-l)\£+l 

2X£*+(m-£)X£ 

*X-1 


N, 


Let  the  reduced  transfer  function  obtained  in  (4.5)  be  parameterized  as 
follows: 


d  ■  sn  +  bjS0"1  +  b2sn"2  +  ...  +  bn 
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— 1 
<T 
a 
+ 
t-* 

_ 1 

bn+2 

bn+n 

b2n+l 

b2n+2 

b3n 

• 

• 

n-1 

• 

n-2  . 

• 

• 

s  + 

• 

s  +  . . .  + 

• 

b4n+l 

b£n+2 

b(i+l)n 

b(4+l)n+l 

• 

• 

- 

b(a+l)n+2 

• 

• 

bU+2)n 

• 

• 

• 

• 

• 

_b(m+l)n_ 

Let  a.  *  (a^.ag . ax*(m+£)^‘  *s  c^ear  from  (4-5)  that  each  bj  is  a 

function  of  a^  and  will  be  indicated  as  such  b.(aj  whenever  necessary. 

J 

It  is  in  fact  a  polynomial  in  a.. 

Let  f:  R"^*)  — ►  Rn<m+1)  be  defined  as: 

'  (bl’b2 . Vl)n>  *  * 

The  function  f  therefore  describes  the  “transformation"  precisely.  The 
structure  of  f  win  play  a  very  Important  role  In  the  proof  of  the  main 
result  where  It  will  be  required  to  show  that  the  Jacobian  of  f,  Is  ' 
full  rank  at  some  point  a..  The  following  Proposition  addresses  this  issue. 

Proposition.  Let  J^  be  the  Jacobian  of  f.  There  exists  a  point  a.  (i.e. 
a  specific  transfer  function  of  the  type  (4.1))such  that  J^(a)  Is  full 
rank. 

Proof: 

The  minors  of  $  were  Introduced  in  (4.4). 


Oil  "  9 


T  Mil  3 


T  i.i  T  || 


+  ...  + 


1 » (t-l)x 
h2,(i-l)x 


h  .(l-l)X-l  .  .  h 

hils  •••  hi,(i-l)x  * 


Let  A  be  the  nxn  matrix 


"  1  0  0 

0  0  ...  0 

■»n  ha  hu 

1  0  ...  0 

**i2  "a  "a 

•  •  * 

hn  hn  hu 

•  •  0 

•  •  • 

*MM)x  •'a.ti-Di  h».(t-U» 

•  •  « 

•  •• 

•  •• 

0  0  0 

hl.(»-l)»  *2,(»-l)x  h*U 

m 

0  0  0 

•  o  • 

0  0  0 

•  •  • 

•  o  • 

0  0  0 

•  •  • 

0  0  0 

1  o 

•’ll  ha 


SiU-m  "a.U-m 


&  block  coIwm. 


Let  <xj 


(al»*2»  •••  ^  Dj^p 

Wi+l*  •'*  a2XJ^  2°d  ro*  of  DRP 


%+i 


ln(mn-l)+l  an(m+i))  last  row  of  NRP* 
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Let 


; 


,  an  U+l)n  x  (t-l)n  matrix 


,  an  (m-£)n  x  (m-i)n  block  diagonal  matrix. 


Let  B  be  an  (a+l)n  x  (fc-l)n  matrix  and  E  some  (m-«.)n  x  (i-l)n  matrix. 
As  will  be  seen  very  shortly  the  structure  of  B  and  E  is  not  needed 
for  the  proof. 

Then  is  given  by 


which  can  be  transformed  by  a  similarity  transformation  to: 


q+*|  y  *  **  iuii  i  ailN  DIIU  wile  vvi  iea|Aiin«iii^  n\«/ 

coprime. 


9  ■  \  •»  /  w 


with  a  block  diagonal  (square) matrix  with  A's  on  the  diagonal. 
It  is  easily  seen  that  for  the  point  a.  defined  by: 


^  -1  0 

0  0-1 


D  «  IsA  ,  N  -  : 


1  1  1 


A  Is  full  rank  and  so  therefore  is  Jf.  This  completes  the  proof  of 
the  Proposition. 


The  stage  has  now  been  set  for  stating  and  prooving  the  main  result 
of  this  paper. 


Theorem 


Let  P  *  nrpdrp_1  be  an  1,1X4  (m21)  strictly  proper  transfer  function 


where 


drp“IsX  +  da-isX_1  +  -+do 
NRP  ■  Px-f  1  +  1"Nol* 

M  N 

Let  m|n  and  qiO  a  fixed  Integer.  Let  k  *  (q+l)m  +  qsn  +  q*1,  *(s)  the 


closed  loop  characteristic  polynomial  and  define: 

«VDRP>tRln*‘)n|°IW>  -  Is*  *  D^jS*-1  .....  00,  NRp  -pij.f*-1  *  ...  « 

k,  -  Ki  h 

S  ■  {(sj,  s2»  ....  sk)CR s1  real} 

JfN  n  si  „  p(m+*)n  v  ok  For  which  there  exists  a  proper  compensator  of 

l\r.RP,uRp,^  ck  x  k  or<jer  q  such  that  Si|S2>>j>fSk  are  roots  of 

Then  Z  Is  a  generic  subset  of  R(m+*)n  x  Rk. 
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Remark:  The  requirement  that  m|n  is  introduced  merely  for  convenience. 
Proof: 

From  (4.6)  let  y.  =  [Yq,  xq-1,  ...,  xQ,  YQ].  Then  as  in  the  mxl  case 

1  •  Sq+1(d,N)Q  =  -Csj’.Sg1 . sj]  -  a  •  Q  .  (4.7) 

This  follows  directly  from  the  reformulation  of  the  mx£  problem  as  an  mxl 
problem.  ^q+1(d,N)  is  a  matrix  whose  entries  are  polynomials  in  a_.  Now 
the  set  EC  R^m+a^n  x  Rk  for  which  a  solution  to  4.7  exists  and  is  such 
that  the  corresponding  x(s),  y(s)  are  left  coprime  is  a  Zariski  open  set. 

For  it  to  be  generic  it  must  be  shown  to  be  non-empty.  The  Main  Lemma 
guarantees  this  to  be  true  for  almost  all  b  x  s.  Since  there  exists  an  a^ 
(by  the  Proposition)  for  which  Jf(a.)  is  full  rank,  then  there  exists  an 
open  set  U  such  that  f(a)  e  U  and  f(R^m+£^n)  3  u,  (by  the  inverse  function 
theorem  Luenberger  1969).  This  means  that  E  contains  at  least  one  point. 
This  completes  the  proof  of  the  Theorem. 

The  Theorem  suggests  that  for  "almost  all"  mx£  transfer  functions  of 
McMillan  degree  n  (and  equal  controllability  indices  A  *  y  )  and  for  "almost 
all"  .s,  (k  *  (q+l)m  +  q)  there  exists  a  proper  compensator  of  order  q  such 
that  Sj.Sg*... »sk  are  k  roots  of  the  closed  loop  characteristic  polynomial 
4>(s).  As  in  the  mxl  case  some  roots  may  be  left  unassigned.  In  the  multi¬ 
input  case  the  possibility  does  exist  for  introducing  additional  parameters 
in  the  compensator  that  can  be  used  to  control  the  remaining  roots. 

One  such  possibility  is  to  modify  the  original  compensator  structure. 
Let  C  «  X-1Y, 
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L  , 

1 

! 


J 

» 


x(s) *  0,  ... ,  0 

”y1(s),  y2(s),  ...,  y£(s),  ....  ym(s)“ 

0 

0 

N 

X 

h-1 

Y  * 

*!  L  0 

r 

0 

o 

where  L  is  an  £-1  x  4-1  constant  matrix  containing  (4-1)  free  parameters 

comprising  a  vector  c.  The  compensator  C  ■  X_1Y  is  still  proper  of  order  q. 

Following  a  similar  proof  to  the  one  given  one  can  show  that  k  poles  can  be 

2 

arbitrarily  assigned  for  almost  all  c  x  a  x  s  e  R^'1^  x  x  Rk. 

This  means  that  for  "almost  all"  choices  of  a^  and  _s  a  proper  compensator 
of  the  form  given  above,  which  is  parameterized  by  c ,  (valid  for  "almost  all" 
c),  can  be  constructed  which  assigns  k  poles  of  the  closed  loop  system  to 
s.  *  (s^Sg,. . .  .sk) .  The  freedom  afforded  by  the  presence  of  these  parameters 
can  then  be  used  to  "control"  the  location  of  the  remaining  unassigned  poles. 
An  illustrative  example  is  given  in  section  5. 

Remark:  As  in  the  mxl  case  the  number  of  assignable  poles  increases  as  a 
function  of  q  in  such  a  way  that  if  q  *  y-1  then  all  the  closed  loop  poles 
are  arbitrarily  assigned.  Brasch  and  Pearson  1970  show,  in  an  entirely 
different  way,  that  for  a  controllable  observable  system  adding  a  p-1  order 
compensator  is  sufficient  to  ensure  arbitrary  pole  assignment. 
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5.  Complex  Poles 

The  concern  thus  far  has  been  the  arbitrary  assignment  of  a  number  of 
real  poles.  The  results  remain  valid  for  the  case  when  real  and  complex 
conjugate  roots  are  desired.  It  is  evident  that  generically  a  solution  to 
(3.1)  will  exist.  The  only  requirement  is  the  invertibility  of  Sq+1  •  Q. 
Since  the  solution  contains  the  compensator  parameters  the  critical  issue 
is  to  show  that  the  unique  solution  is  real .  This  issue  is  addressed  in 
the  following  Lenina. 


lemma. 

Let  i  *  (S1*S2*S3»S4>. • • *s2j-l*s2j,s2j+l . sk> 

be  a  set  of  k-2j  real  and  2j  complex  conjugate  values,  ((Sj^) »  (s3,s4)  ... 

( s 2 j _ i » s 2 j )  are  J  complex  conjugate  pairs  and  s2j+i»****sk  k_2j  red  va,  iS- 
The  unique  solution  ^  of  (3.2) 

£  •  Vl(d,N)  '  Q  ‘  '  . . k]  '  ’  Q 

(whenever  it  exists)  is  real. 


Proof; 

Let  the  1th  row  of  be  thought  of  as  the  coefficients  of  polynomial 
fcj .  Then 

♦l ( Si ) ( Sg) » 


Vi  • « 


♦2^si^ »^2^s2^ *  ****  *2^sk^ 


L*k(*i,**k(s2,»  ****  *k^sk^ J 


M 
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4 
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Since 

Is  -  [sj.  sj,  ....  sJjM"1  -  a  •  QM"1 
then  ^will  be  real  if 

[sf,  Sg.  •••.  -.[u1#  u2»  ....  uk]  is  real  for 

every  integer  t  2  0. 

If  M..  is  the  i,j  minor  of  M  then 

•  J 

u  [(-D^sX-1  +  (-DJt2^|2  *  -  * 

J  (-l)J+^j(s1)Mj1  ♦  t-l)J+5*j(s2)Mj2  +  ...  ♦  <-l)J+k*j(sk)Mjk 

where  det  M  is  expanded  using  the  jtfl  row. 

It  is  not  difficult  to  see  that  Uj  *  u^*  where  *  indicates  complex 
conjugate.  This  means  that  Uj  and  therefore  £  is  real. 

Using  this  result  one  can  easily  see  that  the  Main  Lemma  and  Theorem 
still  hold  if  js  contains  real  and  complex  conjugate  values. 

The  following  example  helps  to  illustrate  the  pole  assignment  method 
in  the  multi-input  multi-output  case. 


Example  2 

Let  m-2,  £*2,  n*4  and 


RP 


f1  o“ 

? 

"o  oi  ro  -r 

m 

1 

s  nrp  " 

s  + 

1 

rH 

o 

Kr 

_i  oj  L1  i- 

Using  the  modified  compensator  structure  given  In  (4.8), 

Y 


'1  O’ 

1. 


"*1  H 
_0  c  J 


this  transfer  function  is  reduced  to  the  2x1  transfer  function 


d 


s  +  cs 


"0“ 

0 

r  °i 

0 

■c“ 

~  c~ 

N  * 

_1_ 

s3  + 

LiJ 

sz  + 

_0_ 

s  + 

.0. 

The  results  suggest  that  with  a  constant  compensator  (q*0),  2  poles  can 
be  assigned. 

Let  s2  =  (-1  +  j2)  ,  s2  *  (-1  -  j2)  . 


The  compensator  parameters  y^,  y2  are  given  as  the  solution  of 


[yj^  s 

0  0  c  c 

r  3  3i 

S1  s2 

a  -  [s14,s24]  -  [0,c,0,0] 

r  3  3-1 

S1  s2 

i  1  0  0 

2  2 
S1  s2 

2  2 
S1  s2 

S1  s2 

S1  s2 

1  1  _ 

i  1 

Computing  the  solution  yields 


*1 


(c-3)25 
8c  * 


9 


and  the  compensator  in  parametric  form 


(c-3)25 

8c 


7+3c 

8 


This  compensator  makes  the  closed  loop  characteristic  polynomial  equal  to 

♦(s)  «  (s2  +  2s  +  5) (s2  +  |(c-3)s  +  |(c-3)). 

One  can  easily  see  that  for  "almost  all"  choices  of  c(i.e.  c^O)  the  compen¬ 
sator  makes  -1  +  j2  and  -1  -  j2  two  of  the  closed  loop  poles.  In  this 
simple  example  the  remaining  two  roots  can  be  explicitly  expressed  as 
functions  of  c. 


In  particular  for  c  >  3,  S3  and  s^  are  guaranteed  to  be  stable.  Had  one 
used  the  compensator  structure  suggested  in  (4.3),  c*l,  it  would  correspond 
to  s3  *  1.5542476,  s^  *  -.8042476.  which  includes  an  undesirable  pole. 

The  above  suggestion  becomes  a  very  powerful  tool  in  that  the  compen¬ 
sator  C  is  given,  parameterized  by  _c,  that  assigns  k  of  the  closed  loop 
poles.  Since  the  remaining  poles  in  general  depend  on  ct  they  can  in  turn 
be  controlled. 
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6.  Conclusions 

Using  an  approach  involving  Sylvester  Resultants  it  is  demonstrated 
that  generically  min(n+q,  )nH-q)  closed  loop  poles  can  be  arbitrarily 
assigned  with  an  output  fedback  compensator  of  order  q.  It  Is  further 
suggested  how  the  locations  of  the  remaining  unassigned  poles  could  be 
controlled.  The  approach  is  different  than  the  ones  followed  by  Antsaklis 
and  Wolovich  1977,  Brasch  and  Pearson  1970,  Kimura  1975,  Davison  and  Wang 
1975.  For  the  appropriate  cases  the  result  is  an  improvement  of  the 
earlier  result  CAntsaklis  and  Wolovich  1977,  Kimura  1975)  for  dynamic 
output  feedback.  The  method  of  solution  can  be  easily  programmed  on  a 
digital  computer. 

It  is  my  belief  that  the  results  in  the  multivariable  case  can  be 
strengthened  by  exploiting  more  effectively  the  compensator  structure. 
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